Orbital-spin order and the origin of structural distortion in MgTi 2 04 
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We analyze electronic, magnetic, and structural properties of the spinel compound MgT^CU using 
the local density approximation+U method. We show how MgTi2C>4 undergoes to a canted orbital- 
spin ordered state, where charge, spin and orbital degrees of freedom are frozen in a geometrically 
frustrated network by electron interactions. In our picture orbital order stabilize the magnetic 
ground state and controls the degree of structural distortions. The latter is dynamically derived 
from the cubic structure in the correlated LDA+U potential. Our ground-state theory provides 
a consistent picture for the dimerized phase of MgTi2C>4, and might be applicable to frustrated 
materials in general. 

PACS numbers: 71.30.+h, 72.80.Ga, 61.66.Fn, 71. 15. Mb 



I. INTRODUCTION 

Over the last years, a lot of progress has been achieved 
in the study of transition metal oxides on frustrated lat- 
tices. Interest in these systems stems from the richness 
of their novel properties: the unexpected variety of or- 
dered states and transitions between them, and the com- 
plexity of the underlying physics. Transition metal oxide 
(TMO) spinels AB2O4, with magnetic B ions forming a 
pyrochlore lattice, give the unique possibility to explore 
how the natural tendency of correlated systems to de- 
velop magnetic, orbital, and charge order is effected by 
geometrical frustration^ 2 - The best studied example of 
TMO spinels, and historically the first one, is magnetite 
Fe304 which shows a high Curie temperature and under- 
goes the Verwey transition at Ty ~ 120 K<2 

Recently, B-spincl MgTia04, which is characterized by 
a pyrochlore lattice of Ti 3+ magnetic ions with one sin- 
gle electron in the ^g-uianifold, has attracted much at- 
tention due to a very peculiar phase transition from a 
metallic to a spin-singlet insulating phase near 

Tc ~ 260 K. 3 The signature of the insulating state is 
the optical gap (« 0.25 eV at T=10 K) observed in optical 
conductivity spectra** Taken together, dc-resistivity 3,5 
and optical^ measurements consistently indicate that 
MgTi204 undergoes a sharp metal- insulator transition 
(MIT) with no sign of Drudc weight at low frequencies 
below Tc- The MIT is accompanied by a structural tran- 
sition from cubic to tetragonal symmetry, with a con- 
comitant drop of the magnetic susceptibility and a re- 
sistivity jump below Tc- Neutron diffraction and X-ray 
measurements 7 - indicate spin-dimerization: Ti-Ti dimers 
are formed in a helical pattern in the spin-singlet state. 
Further structural refinements^ reveal that Ti-ions move 
away from the center of the TiOe octahedron. In the 
low-T phase two out of six Ti-Ti bonds gets closer, sug- 
gesting the formation of chemical dimers. These findings 
have suggested a removal of the pyrochlore degeneracy 
by a ID helical dimerization of the spin pattern, with 
spin-singlets (dimers) located at short bonds. 




FIG. 1: Orbital order in the tetragonal crystal structure of 
MgTi 2 4 within the LDA+U framework ({7 = 3 eV). The 
electron density of the most populated t2 g orbital is plotted. 
Lattice constants are taken from Ref.— . The orange bonds 
represent the shortest bonds in the distorted structure, and 
the corresponding atoms are displayed in red (Ti), blue (O) 
and white (Mg). 



There were several theoretical attempts to understand 
the nature of spin-dimerization and the origin of the MIT 
in MgTi204. Khomskii and Mizokawa 9 assumed that the 
system is close to an itinerant state, and explained the 
formation of the nonmagnetic spin-singlet state by ex- 
ploiting the concept of a orbitally driven Peierls state, 
leading to the formation of quasi ID bands. In this 
case, the magnetic changes across the structural tran- 
sition in MgTi204 can be understood within the picture 
of ID Peierls transition driven by the ordering (Fig. [T|) of 
the d xz and d yz orbitals. This picture is consistent with 
B3LYP functional (GGA) calculations 7 for the tetrago- 
nal phase of MgTi204, showing that the xz and yz or- 



bitals are occupied while the xy states are pushed to high 
energies. 

Another approach is based on the assumption that the 
low-tempcraturc tetragonal phase of MgTi2 04 is Mott in- 
sulator. 10 In this case, the ground state of MgT^C^ can 
be found by studying an effective low-energy spin-orbital 
super-exchange Hamiltonian. It was shown iniS that the 
orbital degrees of freedom in MgTi204 modulates the 
spin exchange couplings, providing an explanation for the 
helical spin-singlet pattern observed in»i It appears that 
the minimum energy configuration corresponds to such 
an orbital ordering for which the maximum number of 
spin-singlet dimers is formed. However, there are vari- 
ous dimer coverings of the pyrochlore lattice which have 
the same energy. As in spin-Peierls systems, the increase 
of magnetic energy gain can be achieved by the shorten- 
ing of bonds, where dimers are situated. Therefore, each 
type of dimer coverings corresponding to a particular or- 
bital ordering induces a different distortion of the lattice, 
which costs a different elastic energy. The ground state is 
chosen simply by finding a minimum energy state, which 
in case of MgTi2C>4 corresponds to the minimal enlarge- 
ment of the unit cell. 

However, no theoretical studies have been yet per- 
formed in order to understand the origin of the MIT and 
the formation of the spin-singlet state. Is it the intrinsic 
lower dimensionality of MgTi2C>4 that causes the forma- 
tion of the spin-singlet state, or is the driving mechanism 
of different type? 

In this work, in order to address these questions, we 
perform quantitative investigation of the band structure 
of MgTi2C>4. Using local (spin) density approximation 
plus Hubbard U, L(S)DA+U, approach^ we show how 
an explicit incorporation of electronic correlations allows 
for a realistic description of the insulating ground state of 
MgTi2C>4. Our results reveal that electronic correlations 
are fundamental in stabilizing the dimerized ground state 
of MgTi2C>4. We show that MIT is driven by correlation 
(via LDA+U) induced orbital order (00) rearrangement. 
Therein, U controls band splitting towards an orbital- 
insulating state without full orbital polarization (OP). In 
our picture 00 stabilizes the spin-singlet ground state, 
which in turn controls the degree of structural distor- 
tions. This finding is consistent with the superexchange 
spin-orbital description^ The opening of the electronic 
band gap is understood in terms of an orbital-selective 
MIT 12 on a quasi ID network. 



II. RESULTS AND DISCUSSION 

To elucidate the interplay between spin, orbital, and 
charge degrees of freedom in MgTi204 we perform cor- 
related scalar-relativistic band-structure calculations us- 
ing local (spin) density approximation plus Hubbard U, 
L(S)DA+U, approach^i We employ the linear muffin-tin 
orbitals (LMTO) scheme in the atomic sphere approxi- 
mation, with combined correction termsJ^ 
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FIG. 2: LDA+U orbital resolved DOS for the tetragonal 
(P4i2i2)£ phase of MgTi 2 4 : U = eV (top left), U = 3 eV 
(top right), (7 = 5 eV (botton left), and U = 6 eV (botton 
right). Notice the orbital selective nature of the correlation 
induced metal-insulator transition between 5 eV < U < 6 eV 
solutions. 



Self-consistency is reached by performing calculations 
on a 12x12x10 k-mesh for the Brillouin zone integra- 
tion. Experimental atomic positions^ for the tetragonal 
(P4i2i2) low-T phase were used in the calculations. The 
radii of the atomic spheres were chosen as r = 2.4789 
(Mg), r = 2.0022,2.0106 (01, 02), and r = 2.5857 (Ti) 
a.u. in order to minimize their overlap. 

Our results for the paramagnetic (LDA+C/) phase of 
MgTi 2 04 are presented in Fig. [2j Therein the evolution 
of the 3d 1 correlated density of states (DOS) is shown. 
In cubic spinel oxides the <i-electron orbital sector splits 
into the low-energy ti g and high-energy e g orbitals. Be- 
low we will denote the ti g orbitals in the tetragonal phase 
as XY, XZ, and YZ, which differ from the d xy , d yz and 
d zx in the cubic phase (see details ini£) . The correspond- 
ing DOS are plotted in red (XY), orange (XZ), and blue 
(YZ) in Fig. [2] The DOS with dominant e g contribution 
are shown by green and cyan curves. Due to partial oxy- 
gen relaxation, the crystal field has also a small trigonal 
component, which further splits the ti g manifold into a\ g 
and e' g sectors .— We note that the XZ orbitals directed 
along one of the short Ti-Ti bonds can only be formed as 
a linear combination of the a± g and e' g orbitals. Accord- 
ing to LDA results for cubic MgTi204 the occupation 
of the e' g states is appreciably higher than that of the 
ai g one. Already in LDA calculations for the dimerized 
phase, the degeneracy of the e' g orbitals is lifted as a con- 
sequence of the local distortion of TiOe octahedra and 
overall tetragonal symmetry One of them (orange curve 
in Fig. 2) acquires XZ character and becomes more pop- 
ulated than the Y Z one (blue). The XY orbital (red), 
originating from the a\ g state, forms a peak just above 
the Fermi level (Ep), its occupation being less than that 
of the other two. This picture does not change much in 
the weak correlation regime (U <3 eV). 

With increasing U (3 eV < U < 5 eV) transfer of 
spectral weight significantly modifies the weakly corre- 
lated scenario such that the XZ orbital is now almost 
half filled (nxz > 0.8), whereas Y Z and XY states are 
pushed towards the e g bands, reducing the t2 g — e g charge 
transfer gap. In this large U regime the overall electronic 
bandwidth is enhanced, accompanied by a reduction of 
the DOS near Ep and selection of a single orbital chan- 
nel to form the insulating state. We observe a clear ten- 
dency towards a pseudo-gap formation in the XZ and 
YZ channels (U = 5 eV), while the XY orbitals are 
pushed to higher energies. This indicates OP as the pre- 
cursor of the MIT. By further increasing U we obtain an 
insulating state at U — 6 eV. Our paramagnetic insu- 
lating solution is characterized by an almost half filled 
XZ band, with appreciable splitting between occupied 
bonding and unoccupied antibonding states. 

The bonding-antibonding splitting observed at large 
U > 5 eV (Fig. aids the formation of the charge gap. 
The insulating state can be viewed as a consequence of 
strong hopping between the XZ orbitals along the short 
Ti-Ti bonds, leading to robust singlet character of these 
bonds (see discussion below) along the c-axis. Thus, ac- 
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FIG. 3: Orbital ordering of MgTi2 04 as obtained within 
LSDA+U (U=3 eV) in the antiferro-pseudocubic (a, upper 
panel) and antiferro-dimerized phase (b, lower panel). Red 
and blue colors denote the ground-state orbital with different 
(t, 1) spins. 



cording to our results MgTi 2 04 undergoes an orbital se- 
lective MITfi^ which is caused by correlation assisted or- 
bital rehybridization»i^ 





FIG. 4: LSDA+U orbital resolved DOS for U = 3 eV in 
the antiferro-pseudocubic phase (upper panel) and in the 
antiferro-dimerized phase (lower panel) of MgT^CU. Within 
dimers, DOS of different sites are mirror images with respect 
to spin space. 



The changes in 00 and crystal structure are shown in 
Fig. [T] Consistent with- ferro-orbital order (FOO) along 
the short Ti-Ti bonds, which stabilizes the molecular 
bond (orange) formation, starts to develop already for 
U < 3 eV. The spiralling of the dimer bonds along the 
c-axis is visible in Fig. [TJ For U — 3 eV we obtain FOO, 
with the ground state orbital pointing through bond's 
direction. 

We now turn our attention to orbital, spin and charge 
responses in the magnetically ordered state. In Figs. [3] 
and |4] we display our LSDA+U ({7 = 3 eV) results for an 
idealized pseudocubic (upper panel) pyrochlore structure 
and the low-T tetragonal structure (lower panel) Z 

In the idealized pseudocubic (PC) phase all atoms are 
arranged as in the high-T cubic structure supplemented 
by a reduced (P4i2i2) tetragonal space group. This com- 
bination of atomic arrangements and crystal symmetries 
allows for the discovering of novel magnetic and orbital 
reorientations in the undistorted pyrochlore structure of 
MgTi204, which are in principle allowed in the vicinity 
of the MIT points 



FIG. 5: Additional magnetically ordered LDA+U solution 
({7 = 3 eV) with the same FOO as the ground-state solution, 
however with FM spin order in the a, b plane and AFM spin 
order along c. Red and blue colors denote the orbitals with 
different (f , J.) spins. 

The orbital resolved DOS of the PC structure (Fig. @J 
shows large OP and a small charge gap between the 
ground (XZ) and the first excited XY orbital. Notice 
the dramatic rearrangement of the minority states of the 
ground state orbital, which are shifted to energies above 
3 eV. 

Within the tetragonal metric of the low-T phase we 
find a similar evolution for the magnetically-ordered elec- 
tronic DOS as for the pseudocubic regime. 

In the lower panel of Fig. [3] we display the ground-state 
00 in the distorted phase of MgTi 2 4 . The LSDA+U 
solution gives the same antiferromagnetic arrangement of 
spins along the short Ti-Ti bonds and in the a, b plane of 
the pyrochlore lattice, as found for the PC structure. The 
lower panel of Fig. 3 shows, however, that the canting of 
the XZ orbitals away from the Ti-Ti bond plane becomes 
signicantly smaller when the low-T structural distortions 
are taken into account. Also, the average Ti-0 distance 
in the XZ plane (2.068 A) is appreciably larger than in 
the YZ (2.052 A) and XY (2.041 A) planes, which means 
that the occupation of the XZ state becomes energeti- 
cally more favorable. This suggests a novel scenario for 
quasi ID chains^ in the strongly frustrated network of 
MgTi204, where the crystal-structure itself gets modi- 
fied by the onset of 00. The small changes in the atomic 
positions of the Ti-ions shown in Fig. [3] additionally sug- 
gest that crystal structure transformations are coupled 
to the correlation-induced orbital re-orientation. 
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FIG. 6: ([7 = 3 eV). Ti-Ti shortest dimers are represented 
in mauve. Tetramers and spirals ( only one representative is 
shown in orange and red, respectively) arises from distorting 
the pyrochlore lattice (black transparent bonds) from cubic 
to tetragonal. Only Ti atoms are displayed for clarity. 

Another magnetically ordered LSDA+U solution 
(Fig. [5]) with the same FOO but FM spin order in the 
a, b plane and AFM one along the c direction can also be 
obtained for both the PC and dimerized structures. The 
order along the short Ti-Ti bonds remains FOO with 
AFM spin order, however, the Ti spins at the ends of 
intermediate Ti-Ti bonds with AFOO (XZ-YZ) in the 
c direction are antiparallel. As a consequence, the total 
energy of this solution is higher than the ground state 
solution shown in Fig. [3l in which the order along inter- 
mediate bonds is AFOO with FM spin order. 

The dimcrization of AFM spin chains and the forma- 
tion of spin-singlets located on short bonds explains the 
crystal structure of MgTi204 , which can be seen as a 
collection of helices running along the c axis. This he- 
lix is formed by an alternation of short-long-short-long 
bonds as found in^ As a consequence of the selective 
bond shortening superstructure features like tetramers 
can be recognized (Fig. [5]). 

The static calculations have provided evidence on the 
role of electronic correlation in the formation of the spin- 
singlet state. To further unravel the coupling between 
structural rearrangement and correlation-assisted 00, 
we have performed full structure relaxation in the corre- 
lated LDA+U potential. Starting from the cubic struc- 
ture with pyrochlore sublattice, Parrincllo-Rahman 19 



(PR) structure relaxation^ was performed. The La- 
grangian formulation of PR allows for constant-pressure 
relaxation simulations, as variation of size and shape of 
the simulation box are allowed. This is typically used 
in connection with pressure- induced phase transitions. 
Here U is acting as some sort of "pressure" on the or- 
bital channels, whose rehybridization brings structural 
changes about. To properly capture this effect it is thus 
important to dispose of a molecular dynamics setup al- 
lowing for unbiased geometry changes. For U = 3 eV 
structural lock-in into the tetragonal dimerized phase is 
achieved with selective shortening of a subset of bonds 
{molecular dimers) as displayed in Fig. [6] The role of lo- 
cal Coulomb interaction U as symmetry-reducing agent 
through selection of one orbital channel is apparent in 
the frame of the dynamic calculations. The found di- 
mensional reduction selects thus a particular topological 
order among the various " degenerated valence bond con- 
figurations" , allowed in the quantum dimer framework of 
frustrated lattices^ in agreement with experiments. 



III. CONCLUSION 

In conclusion, we have studied the ground state or- 
bital, charge and magnetic properties of MgTi 2 04 using 
the LDA+U technique. Our results suggests an orbital 
selective^ picture for the metal-insulator transition of 
MgTi204, which is driven by large spectral weight trans- 
fer due to correlations (via LDA+U) . This in turn intro- 
duces a new type of bonding-antibonding splitting near 
the Fermi energy, which is characterized by large orbital 
polarization of excited orbitals. Using the LSDA+U re- 
sults we have computed the orbital order pattern in both 
high-T pseudocubic and dimerized phase. The insulating 
low-T phase is shown to be driven by electron correla- 
tions in concert with spin and orbital order. The exis- 
tence of short Ti-Ti distances along the c-axis follows as 
a consequence of spin dimerization on bonds (molecular 
dimers) . Finally, consistent with x-ray diffraction exper- 
imentsi the resulting superstructure features (tetramers, 
spirals) are shown to be carried by correlation over orbital 
order, explaining magnetic and structural properties of 
MgTi204 in its low-T insulating phase. 
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